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Abstract — In this paper non-group permutation modulated 
sequences for the Gaussian channel are considered. Without the 
restriction to group codes rather than subsets of group codes, 
arbitrary rates are achievable. The code construction utihzes 
the known optimal group constellations to ensure at least the 
same performance but exploit the Gray code ordering structure 
of multiset permutations as a selection criterion at the decoder. 
The decoder achieves near maximum likelihood performance at 
low computational cost and low additional memory requirements 
at the receiver. 

I. Introduction 

The history of permutation modulation goes back more than 
i forty years where it was first introduced by Slepian in [1] and 
'more generally in the framework of group codes in [2]. The 
appropriate definitions are 

Definition I.l 

Let X denote an unit vector in R" and G a finite group 
(order M — \G\) of orthogonal n-hy-n matrices (or a finite 
orthogonal representation of an abstract group). 

1) The pair {G,x) determines the (M, n) group code with 
initial vector x defined as Gx = \gx \ g € G} 
• 2) If G is a permutation group of degree n then the corre- 
] sponding group code is called permutation code 
■ 3) If G is the full symmetric group Sn of degree n then 
\ the corresponding permutation code is called permutation 
modulation 

I 

Since G consists of orthogonal matrices, \\gx\\ = \\x\\ = 1 for 
all g € G, each codeword can be identified with a point on 
a (n — 1) dimensional sphere. Thus permutation modulation 
generates spherical (or equal energy) codes. Usually there are 
two variants of permutation modulation in the literature, but 
here we are concerned with variant I only as defined above''*. 
It is explicitly specified as follows [1]: 



distinct permutations of the components of the initial vector 
X. Thus the size of the code is given as 



(Mi' 



(mi) 



(1) 



where /ii < • • • < /ife, /i'*) denotes i repetitions of the value 
jjL, and n = mi + • • • + rrifc. The code is given by the set of all 

"'Variant II defined in [1] can be obtained from valiant I by applying all 
possible sign changes in the components of each gx, g £ G, whereas now 
the components of the initial vector x are assumed to be non-negative 



M := 
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ml 



(2) 



mi!---m/c! 

where the notation on the right hand side has been borrowed 
from the multi index notation applied to the vector m = 
(toi, . . . , mfc), whereas \m\ := mi + • • • + mk and ml := 
mil . . . TOfc!. In this context x can be interpreted as a multiset, 
i.e. the set of its components with repetitions and the corre- 
sponding permutations are called multiset permutations. The 
repetitions represented by the vector m have been introduced 
to bound the cardinality M of distinct permutations away from 
nl, and therefore to support a variety of code rates 



1 



R=-log^{M) 



(3) 



for a fixed sequence length n. 

In this setting two prominent questions have been asked 
Ql) Given a group G, how to design the initial vector x, 
such that the resulting group code has maximum min- 
imal distance? (Here it is assumed that each codeword 
has equal probability, such that the distance distribution 
characterizes the maximum likelihood detection) 
Q2) How to decode efficiently? (Note that the size Af in (2) 
is still large in most cases, such that maximum likelihood 
decoding is not practical) 
The first question is known as the initial vector problem 
(IVP) and has been addressed in [1] through numerical search 
only. The further history of the IVP covers in particular the 
following items: In [3] it has been shown that for so-called 
full homogeneous components of group representations the 
IVP can be solved. In the more explicit setting of G — Sn 
and prescribed x of the form (1) with the vector m given 
in advance an optimal solution (i.e. determining the vector 
/i = (/ii, . . . , /ifc)) has been obtained in [4]. An algorithm for 
the general case (arbitrary group G) based on mathematical 
programming has been introduced in [5] and refined in [6] by 
means of generalized geometric programming'''. An explicit 
analytic solution of the IVP for permutation modulation (G = 
Sn) has been presented in a hardly recognized paper [8]: Given 
{n,k) the optimal initial vector (1) is determined such that 

•"This method is described in detail in [7], where the authors also remarked, 
that the algorithm does not guarantee to find the global optimum. 



the signal energy is minimized (in particular the vector m 
is not required to be given in advance as in [4]). Due to 
the Lagrangian approach this method finds local minima, but 
parametrizes all possible solutions. 

For the second question there has been developed an optimal 
solution in [1] for the case of permutation modulation: The 
receiver replaces the first mi smallest components of the re- 
ceived vector by fii, the m2 secondary smallest components by 
/U2, and so on. For later reference let us denote this procedure 
as the function SlepianDetect, which takes a sequence of 
length n as input and outputs the sequence just constructed. 
For an arbitrary group G a decoding algorithm has been 
developed in [9], which refines the function SlepianDetect 
by performing an iterative search afterwards. Moreover, in [10] 
ranking/unranking (resp. demapper/mapper) algorithms for 
multiset permutations with respect to lexicographical order are 
presented, which allow to convert information bit sequences 
to codewords and vice versa. 

Contribution: In this paper the following achievements 
with respect to permutation modulation (variant I) correspond- 
ing to the questions Ql and Q2 are presented: 
Al) The solutions [5], [8] to the initial vector problem provide 
only quantized rates corresponding to the cardinalities of 
subgroups of the permutation group 5„ (or the choice of 
the initial vector respectively), i.e. M is exclusively given 
by (2). Since M is usually still very large, only high rates 
according to (3) can be achieved in this way. In this paper 
a permutation modulation construction scheme allowing 
arbitrary rates while maintaining large minimal distances 
is presented. In particular codes with not too large rates 
in possibly high dimensions n will be constructed. These 
new codes do not possess any group (orbit) structure any 
longer, but can be easily constructed. 
A2) A low complex suboptimal decoder for these codes is 
presented. This decoder, though obtained in a completely 
different manner, shares some properties of the algorithm 
in [9] and also with [11]. One part of the algorithm is the 
adaption of the mapping/demapping functions in [10]. 
As a further motivation for the particular code constructions 
presented here serves a certain application of space-time code 
design described in [12], where space-time block codes are 
transformed into spherical codes. In this situation a target 
(space-time code) rate has been specified together with a 
(possibly large) space-time code block length T. This code is 
then transformed into a spherical code with sequence length 
n> T, thus the spherical code rate is scaled by a factor of ^/n 
and one ends up with a quite small rate in a high dimensional 
sphere which does not necessary fulfill the requirement (3). 

II. Rate adapted code construction 

The construction exploits the results obtained in [8] about 
the structure of the optimal initial vector (1)'^': For prescribed 

■^'Note that the solution is based on some Lagrangian optimization technique 
with discretized constraints. Although not mentioned in [8], the proof given 
exploits the convexity of the Lagrangian functional / to adapt the method to 
the discretized case 



k and the minimal distance held fixed the optimal initial vector 
is a minimizer of the energy functional 

k 

E{m, Ai) = ^ TOzmI (4) 

i=l 

The solution is [8] 

= 5— + -11. 

^ , I = (5) 

where [•] denotes rounding to the nearest integer and A > 
0, 77 < are some kind of Lagrangian parameters which 
parameterize the space of solutions unless (n, R) is fixed 
according to (3). The corresponding permutation modulation 
is SnX, where x := 

Let us now fix some target rate R and a corresponding 
code size N = [2"^] ([•], [-J denote the common rounding 
functions to the next greater resp. smaller integer number). 
Then taking a solution (5) with rate greater than R and 
corresponding size M subject to (2) will be the starting 
configuration for the next step. This step utilizes ranking 
algorithms for multiset permutations. The idea is to find an 
ordering of multiset permutations which reflects the Euclidean 
distance relations between the corresponding code points. This 
is roughly achieved by some Gray code ordering. An algorithm 
which lists all M multiset permutations with respect to this 
ordering can be obtained online from [13]. Having this list at 
hand, select N out of the total of M elements of the list as 
equidistant as possible. To this end define A^o € {0, 1, . . . , N} 
to be the maximum number which satisfies \M/N~\ {N — No — 
l)-hLM/iVjiVo < M-1 and set no := \ M/N]{N - Nq-I). 
With this settings pick the first N — Nq elements in the list 
equidistantly spaced by \M/N~\ and the remaining A^o ones 
starting at position no equidistantly spaced by [M/N\ . The 
result is an injective mapping 

mpGrayEncode : {0, . . . , N - 1} — > {0,...,M-1} (6) 

which parameterizes N out of the M multiset permutations 
with largest possible Gray ordering separation. Let us denote 
the set of N selected multiset permutations by A4V and 
the corresponding (A, n) code by C C 5„f. Clearly the 
minimum distance of C is at least as large as the minimum 
distance of SnX and due to the correspondence between 
Gray ordering of permutations and Euclidean distances of 
codewords, the simple parameterization (6) seems promising 
in order to achieve a large minimal distance for C. 

To allow for a low complex decoding at the receiver we 
need a translation table, which translates lexicographic ranks 
to Gray code ranks. The need arises because although we 
have the algorithm [13] which lists all multiset permuta- 
tions in Gray code order, there is no corresponding ranking 
function available (at least to the knowledge of the author). 
A rank function for (multiset) permutations is a function, 
which assigns to each (multiset) permutation a unique number 
in the range 0, ...n!(M) and establishes an ordering of 



(multiset) permutations. The im'erse mapping is called an 
unrank function. For ordinary permutations there exist rank 
and unrank functions with respect to different ordering criteria, 
including lexicographic and Gray code ordering. For multiset 
permutations lexicographic rank and unrank functions have 
been presented in [10, function Demapping/Mapping]''' with 
average complexity proportional to "*:/2. Let us denote the 
lexicographic rank function by LexRank. It establishes the 
required translation table 

LEX2GRAY : {0, . . . , M - 1} — > {0, . . . , M - 1} (7) 

as the inverse mapping of i i — > LEXRANK(7r(i)) when Tr{i) is 
the multiset permutation corresponding to i e {0, . . . , M — 1} 
with respect to the given Gray ordered list and LexRank(7;) 
the lexicographic rank. This mapping can be stored efficiently 
as a Ust of M integer values at the receiver. 

III. Fast near ML decoding 

The transmission of data proceeds as follows. nR infor- 
mation bits are mapped to one of the N codewords. Let us 
assume message i is assigned to codeword x{i) G C given by 
multiset permutation no. i in A4V applied to the initial vector 
X. x{i) will be transmitted through the Gaussian channel, thus 
the receiver gets 

y = y/pnx{i) + w (8) 

where w ~ A/^(0, 1) denotes a white Gaussian noise vector 
and p the SNR at the receiver (since x{i) has unit norm). 
Algorithm 1 presents the decoding procedure in pseudo code. 
The algorithm performs maximum UkeUhood (ML) decoding 



Algorithm 1 DECODE(y) 
1: MLcandidates ^ 

2: for j = 0; j < 2*^"^; j++ do 
3: y{j) ^ Create VARiANT(j) 

4: z{j) ^ SLEPIANDETECT(?/(j)) 
5: MLcandidates <— 

MLcandidates U {NEwCANDiDATE(x;(j))} 
6: end for 

1: i ^ MLDecode{MLcandidates} 
8: return Message no. i 



with substantially reduced number of candidate codewords. 
The careful selection of candidates is the main achievement 
of the algorithm. It utilizes and refines the detection method 
SlepianDetect [1 1 described already in the introduction by 
creating appropriate variants y{i) of the received sequence y. 

Let us go into some detail of Algorithm 1 now and ignore 
the function CreateVariant for the moment. Then the 
loop in line 2 becomes trivial, i = Q, y{Q) = y. Recall 
that the function SlepianDetect would be equivalent to 
ML decoding if we had taken the full codebook SnX for 

'''Note that different from the presentation in [10] all loops and sum 
boundaries have to be decreased by one for an implementation in C, except 
for the upper sum boundary I in function Mapping. Moreover, j has to be 
chosen as the largest I with the property s; > 



transmission. Since our code is a subset, SlepianDetect 
might fail, but the basic idea of Algorithm 1 is that due to the 
Gray like ordering of the multiset permutations the obtained 
candidate is some kind of neighbor (with respect to Euclidean 
distance) of the transmitted sequence (compare the discussion 
in section II). This fact is implicitly contained in the definition 
of the function NewCandidate, see Algorithm 2: It takes the 



Algorithm 2 NewCandidate(2;) 

1: Ir <— LEXRANK(2:) 
2: gr 4- LEX2GRAY(/r) 

3: return [mpGrayEncode {gr)] 

output of SlepianDetect, calculates its lexicographic rank, 

transforms it with the help of Lex2Gray (7) and estimates 
from that number the codeword number by taking the inverse 
of (6), where its domain has been enlarged to the reals (which 
is denoted by here in Algorithm 2'^^). 

Unfortunately the neighborhood assumption is not correct 
in general since transmission errors may occur everywhere in 
the sequence and this is where the function CREATE Variant 
in step 3 of Algorithm 1 enters the stage. Nevertheless Al- 
gorithm 2 provides codeword candidate numbers from scratch 
without going through the list of all codewords in C, thus its 
computational complexity is very low. The final detection step 
7 denotes ML detection with respect to the codewords Usted 
in the set MLcandidates. Since the cardinality of this set is at 
most 2^^'^ the final ML detection has low complexity also. 

At last let us consider the function CreateVariants. 
The idea is, that at least for not too low SNR values, errors 
occur only in a few places. Recall that since each ordered rrii 
received values will be equated to Hi in SlepianDetect, 
thus an error occurs only, if the smallest or largest of them 
is perturbed so badly, that SlepianDetect assigns a wrong 
viilue Hj, i ^ j to it. So this defect interchanges the com- 
ponents nii and + 1 in the sorted sequence. The function 
CreateVariant loops over all fc — 1 such places and inter- 
changes the corresponding components, compare Algorithm 3. 
This strategy obviously approaches ML performance when the 
SNR grows. 



Algorithm 3 CREATEVARlANT(j) 

1: create binary representation j = X]f=i bi2'-~^ 
2: V; 1 6,=i: interchange component no. m; and m; -|- 1 in the 
sorted version of y 



IV. Simulations and discussion 
The following codes have been constructed: 
1) A code with sequence length n = 25 with N = 323 code- 
words (supporting a target rate of i/s) out of M = 600 

'^'in the spirit of [11] the rounding in the last step in Algorithm 2 
corresponds to the decision regions provided by an appropriate lookup table. 
However, the use of the Gray code ordering provides fast access to the location 
of z in the lookup table 



multiset permutations corresponding to m = (1,23,1) 
and^=(-^,0,^). 

2) This code has sequence length n — 50, N = 1024 
(supporting the target rate 1/5) out of Af = 2450 multiset 
permutations corresponding to rn = (1,48,1) and /i = 

(-TI'O'Tf)- 

3) The last code has sequence length n = 100, N = 1024 
(supporting the target rate 1/ 10) out of M = 9900 multiset 
permutations corresponding to to = (1,98,1) and fi — 

(-Tl'O'Tl)- 

Their performance is shown in Fig. 1 for the low SNR 
regime. 
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n^25, m^[1,23,11, N^323, ML detection 
ri-25, in-[1,23,11, H-323, Lex2Gray detection 
n^5Q, m^[1,48,11, N^1024, ML detection 
n^50, m^[1.48,11. N^1024, Lej;2 Gray detection 
n=100, m=[1,98,1|, N=1024, ML detection 
n-wa, m=[1,9S,1|, N=1024, LexJCSray deter- oi 
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Fig. 1. Performance of codes 1-3 with ML detection and detection 
corresponding to Algorithm 1 
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5, ri=|1.23.1], N=323, ML detection 

5, m=[1,23.1j, IJ=a23, Lfc<2Gray<lttefilian 

3,in=(1,48.1], N=1024, MLdetecfcn 

3, ri=(1.4a.1]. N=1024, Lex2Gray delection 
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Fig. 2. Performance of codes 1 and 2 in the mid SNR regime with ML 
detection and detection corresponding to Algorithm 1 

The mid SNR regime for the codes 1 and 2 is shown 
in Fig. 2. The last simulations approve the expectation that 
the proposed decoding algorithm approaches ML performance 
when the SNR grows. In fact the performance curves match 
quite well. 

Observe that ML decoding amounts for code 3 1024 com- 
parisons of 100 dimensional vectors, while Algorithm 1 singles 
out only 4 candidates for ML decoding which are obtained by 
simple element operations. The complexity of Algorithm 2 is 



approximately nk/2 = 150, thus the complexity of the whole 
decoding procedure (Algorithm 1) scales with nlog(n) (due 
to the sorting operations involved) and is independent of N 
in contrast to ML decoding which scales with Nn, whereas 
usually N n holds. As an illustration of computational 
complexity, table I lists the CPU time for the whole simulation 
process of codes 1-3 in the SNR range from -15dB to -2dB. 
Each simulation has been performed on the same machine 
applying Algorithm 1 and ML decoding respectively. The table 
reveals 97.5 up to 99 percent saving. 





code 1 


code 2 


code 3 


Algorithm 1 


.38 


2.22 


128.54 


ML decoding 


14.73 


280.46 


20346.7 


ratio 


.0258 


.0079 


.0063 



TABLE I 

Comparison of CPU simulation times with and without near ML 

DECODING 



In summary, the encoding method described in this paper 
have been shown to support arbitrary rates with minimum 
distance lower bounded by the optimal IVP in the sense of [8] 
with respect to the underlying permutation modulation. Com- 
plementary the presented decoding algorithm achieves near 
ML performance with substantial savings of computation load 
at the receiver and only small additional memory requirements. 
By the way a practical implementation for determining the 
Gray code rank of multiset permutations has been achieved 
with the help of a simple translation table, which provides 
fast access. Note also, that encoding and decoding both are 
simple in structure. 
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